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In  ihe  NW  Atlantic,  the  temporal  mean  circulation  results  in  large  amplitude  surface  topography 
precluding  the  reasonable  use  of  mean  altimetric  surfaces  or  profiles  for  precise  estimation  of  the 
geoid.  We  examine  a  procedure  which  (II  uses  airborne  expendable  bathythermograph  (AXBTl 
sections  collected  along  altimeter  ground  tracks  to  estimate  the  single-pass  or  instantaneous  surface 
dynamic  topography.  (2)  emplov*  "'ncurrent  or  simultaneous  satellite  altimeter  overflights  to  prov  ide 
a  measure  of  instantaneous  sea  level,  and  (5)  arrives  at  a  precise  estimate  of  the  alongtrack  geoid 
profile  by  differencing  the  instantaneous  altimetric  sea  level  and  the  AXBT-derived  dynamic 
topography.  This  technique  is  applied  to  several  Geosat  Exact  Repeal  Mission  (ERM)  ground  tracks 
in  the  area  of  the  NW  Atlantic  Regional  Energetics  Experiment  (REX)  using  AXBT  survey  data 
collected  in  April  and  July  1987.  Geoid  profile  estimates  are  repeatable  between  these  two  independent 
data  sets  to  within  10—20  cm  rms.  T  his  is  near  the  estimated  noise  level  due  to  barotropic  variability 
in  the  Gulf  Stream  region.  These  AXBT  Geosat-dcrived  profiles  provide  a  more  consistent  estimate  of 
the  alongtrack  geoid  gradient  than  do  available  gravimetric  geoids.  Analysis  of  the  statistics  of  the 
ensemble  of  eollinear  passes  provides  an  estimate  of  the  geographically  uncorrelated  orbit  error  of 
approximately  1-2  m.  The  total  surface  topography  derived  as  the  difference  between  adjusted 
instantaneous  sea  level  profiles  and  the  AXBT  Geosat-dcrived  geoid  profiles  does  not  suffer  significant 
contamination  by  either  geographically  uncorrelated  or  correlated  orbit  error 


I.  Introduction 

The  fundamental  objective  of  the  NW  Atlantic  Regional 
Energetics  Experiment  (REXl  is  a  measurement  of  the  mean 
and  mesoscale  eddy  available  potential  energy  (APEl  and 
kinetic  energy  (KEl  in  the  Gulf  Stream  (from  Cape  Hatleras 
to  beyond  the  New  England  Seamount  C'hain  (NESO). 
Three  emerging  technologies  are  employed  to  meet  this 
objective:  (I)  satellite  altimetry  from  the  U.S.  Navy's  Geo¬ 
sat  satellite,  particularly  from  the  eollinear  Geosat  Exact 
Repeat  Mission  (ERM)  [see  Horn  et  til..  I987|:  (2)  extensive 
airborne  expendable  bathythermograph  (AXBT)  surveys 
and  time  series  from  bottom-moored  arrays  of  inverted  echo 
sounders  with  pressure  gauges  (IES/PGs):  and  (3)  simulation 
of  the  circulation  and  associated  energetics  of  the  Gulf 
Stream  through  primitive-equation  numerical  models  of  the 
N  W  Atlantic.  The  data  collected  as  part  of  REX  arc  designed 
to  be  both  complementary  to  |see  I  hillock  ct  nl..  I989|  and 
easily  assimilated  into  numerical  models  for  ultimate  quan¬ 
tification  of  the  mesoscale  energetics  in  a  dynamically  con¬ 
sistent  framework  (see  Mitchell  et  nl.  1 1987|  for  an  overview 
of  REX). 

Toward  these  objectives  it  is  desirable  to  use  Geosat. 
particularly  Geosat  ERM.  altimetry  as  a  means  of  measuring 
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the  surface  topography  associated  with  mesoscale  structure 
in  the  NW  Atlantic.  Throughout  this  paper  we  adopt  the 
nomenclature  of  Taplev  et  til.  [  1 9S2 1  (see  their  Figure  1)  in 
which  sea  level  relative  to  a  reference  ellipsoid  consists  of 
the  sum  of  surface  topography,  tides,  geoid.  and  errors  (both 
environmental  and  instrumental).  Numerous  investigators 
have  been  remarkably  successful  in  measuring  the  surface 
topographic  variability  associated  with  the  Gulf  Stream 
using  satellite  altimetry  |e.g..  Thompson  et  til..  1983:  //<;/- 
lock  ct  nl..  1989],  However,  the  use  of  the  satellite  altimeter 
as  a  means  of  directly  measuring  surface  topography  (i.e. . 
that  component  of  sea  level,  boih  temporal  mean  and  vari¬ 
able.  due  to  ocean  currents)  remains  elusive  owing  chiefly  to 
relatively  large,  unknown  sea  level  gradients  in  the  back¬ 
ground  geoid  (upon  which  the  topography  is  superimposed) 
and  to  the  limitations  of  the  space-time  sampling  character¬ 
istics  of  the  single  nadir-beam  altimeter.  T  he  assimilation  of 
surface  topography  into  numerical  models  of  the  regional 
circulation  promises  a  solution  to  the  latter  problem. 

Previous  attempts  to  compute  absolute  surface  topogra¬ 
phy  from  satellite  altimetry  have  estimated  the  background 
geoid  using  either  a  gravimetric  surface  \Chcncy.  1 982 1 .  a 
mean  altimeter  surface  or  alongtrack  profiles  {Menard. 
19831.  or  a  hybrid  altimetric  mean-gravimetric  surface  |/.y- 
hanon  and  Croat.  1987).  T  hese  attempts  have  been  less  than 
completely  sntiM.ictory  owing  either  to  the  limited  precision 
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of  the  gravimetric  surfaces  (e.g..  errors  in  the  Marsh-Chang 
gravimetric  surface  are  estimated  to  be  of  the  order  of  the 
Gulf  Stream  surface  topography  [Marsh  and  Chant;.  1 978 J ) 
or  to  the  contamination  of  the  mean  surfaces  by  large- 
amplitude  temporal  mean  surface  topography  [ Marsh  et  at.. 
1984|. 

We  describe  a  straightforward  (somewhat  brute  force) 
process  for  using  AXBT  sections  collected  along  collinear 
altimeter  ground  tracks  in  conjunction  with  concurrent  sat¬ 
ellite  overflights.  This  process  will  provide  alongtraek  geoid 
profile  estimates  which  appear  to  be  very  precise  (repeatable 
to  within  10-20  cm  rms)  and  are  not  contaminated  by  the 
temporal  mean  topography.  Our  primary  use  of  these  geoid 
profiles  is  the  computation  of  surface  topography  from 
Geosat  ERM  sea  level  measurements  along  ground  tracks. 
This  estimated  surface  topography  is  not  absolute  in  the 
strictest  sense  (i.e..  we  are  not  able  to  compute  surface 
topography  relative  to  some  inertial  reference  ellipsoid  tied 
to  the  center  of  mass  of  the  Earth).  Both  long-wavelength 
orbit  errors  in  the  altimetry  and  seasonal  steric  adjustments 
made  to  the  AXBT  data  make  any  attempt  at  such  an 
absolute  estimate  impossible  at  present. 

Section  2  provides  a  step-by-step  description  of  the  pro¬ 
cedure  for  using  AXBT  sections  and  concurrent  (hereinafter 
referred  to  as  ■’simultaneous")  Geosat  ERM  overflights  to 
estimate  alongtraek  geoid  profiles.  Errors  inherent  in  the 
procedure  are  assessed  in  section  3.  In  section  4  the  rate  of 
approach  to  stationarity  by  the  ensemble  mean  sea  level  as 
repeat  revolutions  (revs)  (i.e..  orbits  or  passes)  occur  along 
a  given  ground  track  is  used  to  assess  both  the  order  of 
magnitude  of  orbit  determination  error  (or.  simply,  orbit 
error)  and  its  wavelength  inherent  in  the  Geosat  ERM 
ephemeris.  Conclusions  and  suggestions  for  future  applica¬ 
tion  of  the  procedure  are  discussed  in  section  5. 

2.  ProC  EDURE  FOR  ESTIMATING  GEOID  PROFILES 

A  schematic  diagram  of  the  procedure  used  to  provide 
alongtraek  estimated  geoid  profiles  is  shown  as  Figure  I.  It 
must  be  emphasized  that  herein  we  are  not  directly  inter¬ 
ested  in  improving  the  regional  geoid  of  the  NW  Atlantic 
(i.e..  a  two-dimensional  surface):  rather,  we  are  only  con¬ 
cerned  with  estimates  of  the  geoid  profile  along  Geosat  ERM 
ground  tracks  for  the  purpose  of  computing  regional  surface 
topography  from  Geosat  ERM  sea  level  measurements.  Our 
objectives  are  oceanographic,  not  geophysical. 

Our  original  plan  was  to  collect  AXBT  sections  using  U.S. 
Navy  P-3  aircraft  of  the  U.S.  Naval  Research  Laboratory 
Flight  Support  Detachment  along  roughly  12(H)  km  ground 
track  segments  of  13  Geosat  ERM  ascending  ground  tracks 
over  the  Gulf  Stream  region  of  the  NW  Atlantic  from  Cape 
Hatteras  to  east  of  the  New  England  Seamounts.  We  number 
these  Geosat  ERM  ground  tracks  as  A239  to  A7.  where 
ground  tract  AI  (for  "ascending  I")  is  that  ground  track 
overflying  Bermuda  and  ground  tracks  are  numbered  con¬ 
secutively  to  the  east  along  the  equator.  Three  scpaialc  sets 
of  AXBT  surveys  to  complete  this  purpose  were  carried  out 
(December  1986.  April  1987.  and  July  1987).  Optimistic  plans 
called  for  the  repeat  (in  triplicate)  of  each  of  these  13 
sections.  Problems  with  receivers  on  board  the  aircraft 
prevent  us  from  using  the  data  collccteu  in  December  '986  in 
this  study.  Weather  and  aircraft  problems  also  took  their  toll 
on  the  volume  of  data  collected  in  April  1987  and  July  1987. 
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Fig.  I.  Steps  for  estimating  alongtraek  geoid  profiles  from  si¬ 
multaneous  altimetry  and  AXBT  sections.  Steps  diagrammed  in 
Figures  l<i-l/i  correspond  to  those  plotted  in  Figures  3n-3 h. 


Figure  2  posts  the  actual  AXBT  drop  locations  used  in  this 
study.  Note  that  simultaneous  AXBT  sections  and  Geosat 
ERM  overflights  (as  listed  in  Table  1)  were  obtained  for  nine 
separate  ground  tracks  (three  of  these  with  usable  repeated 
sections).  Descending  (NE  to  SW)  altimeter  tracks  were 
excluded  in  this  study  owing  to  significant  data  dropouts 
over  the  NW  Atlantic  on  these  tracks. 

Before  detailing  the  steps  of  the  geoid  profile  estimation 
procedure  as  diagrammed  in  Figure  I .  and  as  illustrated  by 
example  in  Figuie  3.  we  briefiv  examine  the  pedigree  of  the 
altimeter  input  data  (referred  to  here  as  step  0): 

Step  0:  Production  of  Binned  Sea  Level  From 
fiemiit-FKM  Altimetry  Hot,, 

Altimeter  data  input  into  the  geoid  estimation  process 
consists  of  individual  Geosat  ERM  passes  of  time-binned  sea 
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Fig.  2.  Posting  of  AXBT  drop  locations  during  («)  April  1987 
and  (h)  July  1987.  Sections  are  numbered  by  Geosat  ERM  ground 
track  number.  Al  is  the  ascending  ground  track  that  overflies 
Bermuda. 


level  along  a  given  ground  track.  These  data  are  derived 
from  NORDA  (Naval  Ocean  Research  and  Development 
Activity)  data  records  (NDRs)  of  Geosat  altimetry,  consist¬ 
ing  of  time  tags  and  altimeter-measured  satellite  ranges  (i.e.. 
distance  from  the  altimeter's  antenna  to  the  ocean's  surface) 
and  other  geophysical  and  engineering  parameters,  as  ar¬ 
chived  at  the  Naval  Oceanographic  and  Atmospheric  Re- 


TABLE  I.  Summary  of  AXBT  Sections  and  Revs  Used  in  the 
Estimation  of  Geoid  Profiles  From  Simultaneous  AXBT 
Sections  and  Geosat  ERM  Overflights 


Ground 

Track 

No.  of 
Revs 

Simultaneous 

Rev 

Dale  of  Rev 

Date  of  AXBT 
Section 

A239 

12 

10924 

April  15.  1987 

April  17.  1987 

A240 

17 

10967 

April  18.  1987 

April  20.  1987 

12187 

July  12.  1987 

July  13.  1987 

A24I 

17 

1 1010 

April  21.  1987 

April  22.  1987 

12230 

July  15.  1987 

July  15.  1987 

A242 

16 

17773 

July  18.  198’ 

July  17.  1987 

Al 

18 

luv3b 

April  16.  1987 

April  17.  1987 

A2 

16 

10981 

April  19.  1987 

April  20.  1987 

A3 

’() 

11024 

Anri!  22.  I9S7 

April  22.  1987 

12244 

July  16.  1987 

July  15.  1987 

A4 

19 

12287 

July  19.  1987 

July  17.  1987 

A  5 

17 

10866 

April  II.  1987 

April  It).  1987 

search  Laboratory  (NOARL.  previously  NORDA).  Follow¬ 
ing  a  playback  of  Geosat's  onboard  recorded  data,  the  NDRs 
are  transmitted  in  near-real  time  from  the  Geosat  ground 
station  at  Johns  Hopkins  Applied  Physics  Laboratory  ( APL) 
to  NOARL's  Satellite  Data  Reception  and  Processing  Sys¬ 
tem  (SDRPS;  see  Lybanon  and  Croat  [1987]  for  a  more 
complete  discussion).  As  part  of  the  daily  SDRPS  processing 
of  incoming  Geosat  altimetry,  the  following  steps  are  per¬ 
formed. 

1.  Raw  ranges  are  objectively  edited  for  spurious  values 
(WW.  1986], 

2.  Ten-hertz  range  samples  are  smoothed  to  I -Hz  ranges 
R  using  a  running  mean. 

3.  The  NDR  is  merged  with  an  ephemeris  (i.e.,  files 
consisting  of  satellite  height  r  above  a  reference  ellipsoid  and 
the  latitude  and  longitude  of  the  subspacecraft  point  at  each 
t  min  interval)  using  the  time  tags  in  the  NDR. 

4.  The  1-Hz  ranges  are  subtracted  from  the  satellite 
height  r,  as  interpolated  using  a  10-point  Lagrangian  inter¬ 
polator.  to  provide  a  raw  sea  level  (rjohs)  measurement 
(relative  to  a  specified  reference  ellipsoid): 

r~  R  =  V  ohs 

5.  Open  ocean  tide  heights  tj,  computed  from  the  Scli- 
widerski  [  1980]  tide  model  and  a  wave  height  or  EM-bias  plus 
wave  skewness  correction  totaling  5%  of  the  altimeter- 
measured  significant  wave  height  (SWH)  are  removed  from 
this  raw  sea  level  | Lybanon  and  Croat.  1987]  to  form  a 
synoptic  (hereinafter  referred  to  as  “instantaneous")  sea 
level  along  the  track  In). 

V  =  ’Jobs  ~  Vi  +  0.05  *  SWH 

6.  Linear  interpolation  of  the  alongtrack  instantaneous 
sea  level  17  provides  a  profile  of  sea  level  binned  into  1-s 
(about  7  km  alongtrack)  averaging  bins  (17*). 

These  binned  sea  level  profiles  constitute  the  input  data 
for  the  geoid  estimation  procedure  diagrammed  in  Figure  I . 
A  step-by-step  description  of  this  procedure,  along  with  tacit 
assumptions  made  at  each  step,  is  as  follows. 


Step  I :  Computation  of  Binned  Alonptrack 
Mean  Sea  Level 

An  edited  ensemble  of  binned  sea  level  values  (tj*.  where 
index  i  refers  to  a  particular  pass  or  rev  along  the  specified 
ground  track)  is  assembled  over  the  first  year  of  the  Geosat 
ERM.  using  acceptable  altimeter  passes  from  the  period 
November  9.  1986.  to  November  6,  1987.  Each  of  the 
instantaneous  (i.e..  single  rev)  binned  Geosat  ERM  passes 
provides  a  profile  of  sea  level  whose  value  within  each  bin 
( Tj*)  consists  of 


’ll"-  G  +  Ci  +  <6 

where  G  is  geoid  height.  C,  is  instantaneous  sea  surface 
topography,  and  q,  is  instantaneous  (or  single  rev)  orbit 
error. 

For  the  moment,  all  errors  other  than  instantaneous  orbit 
error  have  been  neglected  (see  section  3  for  an  er'or  analysis 
ol  the  procedure).  Here  acceptable  passes  arc  considered  to 
be  those  that  do  not  have  obvious  problems  in  the  measure¬ 
ment  of  sea  level  (such  as  unexplained  steps  of  several 
meters!).  Of  the  approximately  21  repeats  of  the  ERM's 
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Fig.  3.  Example  of  steps  in  an  estimate  of  a  geoid  profile  from  simultaneous  altimetry  and  an  AXBT  section  for 
Geosat  ERM  ground  track  A3.  In  this  case  the  simultaneous  pass  is  provided  by  rev  12244  (July  16.  1987).  roughly 
concurrent  with  the  AXBT  section  (Figure  3/  collected  on  July  15.  1987).  Figures  iu-ili  refer  to  the  procedure  as 
diagrammed  in  Figure  1.  Though  not  obvious  on  the  scale  of  these  panels,  unadjusted  instantaneous  sea  level  (Figure 
3c)  differs  from  adjusted  instantaneous  sea  level  (Figure  3c)  by  a  bias  of  nearly  3  m  (as  estimated  from  a  linear  least 
squares  fit  to  the  residual  profile  in  Figure  3 </). 


collinear  orbit  during  I  year,  as  few  as  12  individual  revs  and 
as  many  as  21  individual  revs  were  used  in  computing  the 
alongtrack  mean  sea  level  profiles  (see  Table  I).  Obvious 
spikes  (i.e.,  large  excursions  in  apparent  sea  level  which 
occur  over  only  a  few  l-Hz  data  points)  are  removed,  and 
linear  interpolation  is  performed  across  the  resulting  gap  (no 
more  than  a  few  seconds  wide).  True  data  gaps  (i.e..  sea 
level  dropouts  of  more  than  a  few  seconds)  are  eliminated 
from  the  averaging  process  with  care  being  taken  that  no 
resulting  artificial  step  function  is  introduced  into  the  ensem¬ 
ble  mean.  This  editing  is  admittedly  subjective,  but  we 


believe  that  it  provides  cleaner,  more  realistic  mean  sea  level 
profiles  than  do  purely  objective  filtering  approaches,  none 
of  which,  to  dale,  have  provided  a  straightforward,  fail-safe 
means  of  editing  altimeter  sea  level  data. 

The  ensemble  mean  sea  level  H  over  N  repeats  of  a  given 
ground  track  is  computed  in  each  bin  as  the  arithmetic  mean 
of  the  bin: 

v 

H  =  UN  £  (77 ,*>  (I) 

I  I 
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where  rj*  is  the  binned  sea  level  value  for  an  individual 
satellite  pass,  or  rev,  i.  The  ensemble  mean  sea  level  H  so 
formed  is  given  by 

H  =  G  +  Z  +  Q 

where  (» is  geoid  height.  Z  is  mean  topography  (sum  of  (,/N), 
and  Q  is  mean  (or  ensemble  residual)  orbit  error  (sum  of 
q-JN). 

To  this  point  we  have  used  the  terminology  “orbit  error" 
to  refer  only  to  that  component  of  the  total  orbit  error  which 
is  geographically  uncorrelated.  Actually,  total  orbit  error  has 
both  geographically  uncorrelated  and  correlated  compo¬ 
nents.  For  the  moment  we  continue  to  ignore  the  geograph¬ 
ically  correlated  component  and  note  that  from  one  pass  to 
the  next  repeat  along  the  same  track,  which  occurs  244  revs 
later,  instantaneous  estimated  orbit  error  q,  is  uncorrelated. 
Thus  we  assume  that  the  summation  used  to  compute  H  is 
performed  over  enough  uncorrelated  orbits  that  Q  =  0  (i.e. . 
residual  geographically  uncorrelated  mean  orbit  error  will  be 
"very  small":  see  section  4). 

Step  2:  Computation  of  Instantaneous 
Alonptrack  Sea  Level 

The  binned  sea  level  profile  for  that  instantaneous 

pass  occurring  nearly  simultaneously  (/'  =  m:  see  Table  I) 
with  the  AXBT  underflight  survey  is  extracted  from  the 
ensemble.  This  individual  pass  suffers  from  residual  orbit 
error  which  we  estimate  according  to  the  tilt  and  bias  (TB)  of 
this  pass  relative  to  the  mean  sea  level  profile.  We  anticipate 
biases  as  large  as  3  m  in  association  with  the  quick-look 
Naval  Astronautics  Group  (NAG)  ephemeris  used  herein. 

To  reduce  the  orbit  error  in  the  individual  simultaneous 
puss  of  sea  level,  the  residual  ( 17'  m).  where 

Vi  in  ~  V  i '  -  m  ~ 

is  computed  in  each  bin.  Thus 

Vi  m  “  I  L 1  -  m  ~~  Z  I  “  Q ) 

or.  assuming  Q  -  0. 

Vi  hi  —  !£/  -  in  “  Z)  ■+■  q,  .  m 

Since  the  residual  orbit  error  is  very  long  wavelength 
relative  to  the  typical  length  of  an  altimeter  pass  over  the 
Gulf  Stream  region  (i.e..  2000  km),  a  reasonable  estimate 
(see  section  4)  of  residual  orbit  error  for  the  simultaneous 
pass  is  obtained  from  the  linear  trend  (tilt  and  bias)  in  the 
residual  (17')  profile  given  by 

Ilf  in  ~  ^  Vi  in) LONG  “  in  ~~  ^  h.ONCi  4  0//  rn)|.ON(i 

such  that 

Vi'  II,  ~  I  Hi  III  -  <4/  m  ~  ^IsilORT  4  i(ll  mlsHOKT 

"  )(i  ,!,-/■)  SHORT 

Subscripts  LONG  and  SHORT  refer  to  linear  and  higher- 
order  alongtrack  trends,  respectively,  in  the  subscripted 
quantity.  Thus  p'  estimates  only  time-variable,  short- 
wavelength  components  of  the  surface  topography.  Such  an 
approach  is  taken  by  Thompson  et  al.  [I983|  in  which  the 


TABLE  2.  Regression  for  Depths  of  I5°C  and  I7°C  Isotherms 
and  Dynamic  Topography  Based  Upon  CTD  Data  of 
Hogg  et  til.  1 1985] 


T.  °C 

Ai-  m 

D0. 

dyn  cm 

Standard  Error 
in  D.  dyn  cm 

*.  dyn  cm/m 

17 

517 

227 

±7.1 

0.191  ±  0.004 

15 

41 1 

227 

±7.3 

0.156  ±  0.003 

ensemble  mean  sea  level  H  is  used  to  estimate  the  geoid. 
Without  ancillary  in  situ  data,  this  is  often  the  only  option. 

The  ancillary  data  provided  by  the  simultaneous  AXBT 
section  allow  us  to  proceed  in  a  somewhat  different  and  more 
inclusive  fashion.  Here  we  will  use  the  binned  mean  profile 
to  estimate  residual  orbit  error  in  each  of  the  individual  revs 
and  subsequently  perform  a  tilt  and  bias  adjustment  to  each 
of  the  instantaneous  binned  sea  level  (17*)  profiles  along  a 
given  ground  track.  The  best  measure  of  sea  level  along  a 
ground  track  at  the  time  of  the  AXBT  section  is  provided  by 
the  adjusted  simultaneous  pass  of  sea  level  (  tj,*  ,„)  given  by 

V  i  ~  m  ~  Vi-  in  ~  TBj  I,.  —  (I  -  f  :  /„  *  q, .  m 

~  l(i  -  m  ~  ^  I  LONG  -  (</  i  -  mh.ONG 

=  G  +  (£,-*.  (SHORT  +  I  LONG  r  <</i  =  mlsHOKT  <-> 

Thus  the  adjusted  estimate  of  sea  level  along  the  ground 
track  (17 ,!„,)  is  formed  by  the  short-wavelength  instanta¬ 
neous  topography  superimposed  upon  the  sum  of  the  geoid 
and  the  long-wavelength  mean  topography  (ignoring  for  the 
moment  short-wavelength  orbit  error;  see  section  4). 

Step  3:  Estimation  of  Instantaneous  Dynamic 
Topography  From  the  AXBT  Section 

A  complete  description  of  the  December  1986.  April  1987. 
and  July  1987  AXBT  data  and  their  analysis  is  provided  by 
Teaptte  el  al.  [1989],  Dastupue  et  al.  [  1988).  and  Mitchell  et 
al.  f  1 988].  Typical  AXBT  sections  used  in  this  analysis 
represent  tracks  of  648  n.  mi  length  (1200  km)  with  AXBTs 
spaced  evenly  at  12  n.  mi  (22  km)  intervals  along  the  track. 
U.S.  Navy  standard  or  shallow  (400  m)  AXBTs  were 
dropped  over  the  continental  shelf  region  with  alternate  deep 
(800  m)  and  shallow  probes  employed  in  the  deep  water  to 
the  south.  Figure  2  depicts  the  positions  of  AXBT  stations 
from  April  1987  and  July  1987.  respectively. 

The  depths  d  of  selected  AXBT  isotherms  ( 1 5eC.  I7CC) 
were  converted  to  dynamic  height  (3  dbar  over  3000  dbar) 
with  the  equation 

D  =  l),i  *  (<S  -  {5(|l  \ 

where  /)„.  and  y  were  determined  by  linear  regression 
from  historical  (1983)  conductivity-tempcrature-depth 
(CTD)  data  from  the  REX  region  [see  Hop  ft  et  al..  1 985 1 . 
Table  2  summarizes  these  results.  The  AXBT -derived  topog¬ 
raphy  (DI5)  based  upon  the  depth  of  the  1 5°C  isotherm  is 
typically  a  bit  less  noisy  than  that  based  upon  the  I7°C 
isotherm  (i.e..  DI7).  although  the  DI5  profile  may  be  some¬ 
what  shorter  (on  account  of  the  AXBT  s  failure  to  reach  the 
deeper  I5°C  isotherm  in  the  Sargasso).  Generally,  the  two 
profiles  are  very  similar. 

It  is  important  to  note  that  there  is  an  unknown  offset 
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between  the  dynamic  reference  level  i3(HH)  dbar)  and  the 
ellipsoid  to  which  the  altimetric  level  is  referenced.  We 
assume  that  this  offset  is  arbitrary  for  regional-scale  pur¬ 
poses  and  remove  it  by  adjusting  the  AXBT-derived  topo¬ 
graphic  profiles  downward  by  2.3  m  (i.e..  the  mean  DI5 
heights  relative  to  3000  dbar).  In  essence,  we  have  moved 
the  3000-dbar  lev  4  to  an  arbitrary  zero  level  relative  to  the 
altimetry.  By  so  doing,  we  necessarily  limit  the  use  of  our 
estimated  geoid  profiles  to  shorter  (i.e.,  nonlinear  trend) 
wavelength  anomalies  alongtrack.  Thus  surface  topography 
computed  as  the  differenced  between  subsequent  altimeter 
passes  of  sea  level  and  these  estimated  geoid  profiles  will,  in 
a  strict  sense,  represent  only  the  alongtrack  anomalies  in  the 
surface  topography  relative  to  the  local  geoid  profile.  We 
must  be  particularly  careful  that  estimates  of  mean  topogra¬ 
phy  based  upon  these  geoid  profile  estimates  are  reasonable 
(e.g..  as  compared  with  mean  topography  computed  from  a 
circulation  model. 

Step  4:  Estimation  of  the  Binned  Alongtrack 
Geoid  Profile 

The  binned  estimated  geoid  profile  along  each  ground 
track  (</*=„,)  is  computed  as  the  adjusted  instantaneous  sea 
level  ( r;,*  „,)  minus  the  adjusted  AXBT-derived  dynamic 
topography  (10 lg*D^m)  from  the  accompanying  AXBT  sec¬ 
tion: 

9f*m  =  Vi-,,,  -  10 lg*Dj  =  m  (3) 

Having  an  estimated  instantaneous  geoid  profile  (<?*  „,) 
along  each  "calibrated"  ground  track  now  allows  the  use  of 
(2)  to  compute  the  sea  surface  topographic  profile  along  each 
ground  track  for  other  revs  (other  than  /  -•  /«)  as 

4/  s  Vi  -  a?-  in  =  (>  +  (fi')SHORT  +  (^)l.ONG  _  It*  in 
+  (</j)sHORT 
-  1 4 /(.SHORT  +  l^h.ONG 

+  1 4/(1  SHORT  9*--  m  =  G  (4) 

This  surface  topography  represents  the  total  suiface  to¬ 
pography  (barotropic  +  baroclinic)  and  is  estimated  as  the 
sum  of  the  short  wavelengths  of  the  instantaneous  topogra¬ 
phy  C,  and  the  long  wavelengths  of  the  temporal  mean 
topography  /.  along  the  selected  ground  track  (plus  some 
small  residual  orbit  error). 

3.  Error  Anai  ysi  s  or  Geoid  Profile  Estimates 

Random  instantaneous  errors  or  uncertainties  will  appear 
as  errors  in  our  alongtrack  geoid  profile  estimate  rig*  ,„)  as 

rig*  ml  “  {r’(l7,  ml  +  lid  +  E-Vrm»>'-’ 

where  C,  is  the  simultaneous  barotropic  component  of  sea 
surface  topography. 

The  direct  appearance  of  these  uncertainties  is  the  major 
drawback  in  our  instantaneous  approach  as  opposed  to 
estimates  of  the  geoid  profile  based  upon  differences  in 
temporal  mean  sea  level  and  dynamic  topography  which  will 
be  addressed  in  a  later  paper.  In  this  section  we  briefly 
examine  likely  errors  in  both  rj,'  and  /),  as  summarized  in 
Table  3. 


TABLE  3.  Root-Sum-Square  (rss)  Residual  Error  Budget  for 
AXBT-E'stimated  Geoid  Profiles 


Error  Source 

rms  Magnitude,  cm 

Emus  in  Altimetric  Sea  Level 

Range  {r(R,)i 

6-10  | Lebanon  and  Croat.  1987] 

Tides  {r((T),),|} 

1  [Lebanon  andCront.  1987] 

EM  bias  !fI|7j<;vvH),)} 

2  [ Lebanon  ami  Croat.  1987] 

Residual  orbit  error  {(•/, (short) 

8  (see  section  4) 

rss  subtotal  {r(  17,'  )} 

10-13 

Errors  in  Estimating  Total  Topography 

Regression  {r(T>,  _  „,)} 

7-8 

Barotropic  error  {rUj}  ,„)( 

10  [H alloc k  el  at..  1989] 

rss  subtotal  {p(f,  _  ,„)} 

12-13 

Total  rss 

Geoid  profile  error  [rig*  _  ,„)} 

16-18 

Root-sum-square  errors  in  Geosat  ERM  adjusted  sea  level 
<f(tj"))  arise  as 

rig  *  )  =  {f  2(R)  +  p ‘( 17,)  +  f  2( tj swh*  +  q2\' 2 

where 

e(R)  error  in  altimeter  range: 
f(t),)  error  in  Schwiderski  tide  model: 
f(j7swH>  error  in  EM  bias  model: 
q  orbit  error. 

Estimates  of  these  four  errors  are  provided  in  the  following 
discussion. 

Errors  in  satellite  range  ir(R))  arise  from  to  unmodeled 
changes  in  the  index  of  refraction  at  the  13.3  GHz  operating 
frequency  of  Geosat’s  altimeter  due  to  intervening  water 
vapor  and  ionospheric  free  electrons.  During  solar  minimum 
(i.e..  1986-1987),  water  vapor  is  likely  the  largest  single 
source  of  error  in  range,  while  unmodeled  ionospheric 
changes  are  small.  Long-wavelength  (i.e..  defined  as  linear 
over  the  track  segment)  water  vapor  error  in  the  simulta¬ 
neous  pass  will  be  removed  during  the  TB  adjustment,  as 
will  as  any  other  linear  trend  in  the  measured  range.  How¬ 
ever,  shorter-wavelength  (i.e..  higher-order)  water  vapor 
signatures  often  associated  with  the  Gulf  Stream  front  re¬ 
main  unmodeled  and  unobserved  (since  Geosat  lacks  a 
bore-sighted  microwave  radiometer  for  water  vapor  correc¬ 
tion).  From  the  error  budget  of  G.  H.  Born  (unpublished 
report.  1984)  given  by  Lebanon  and  Croat  |I987|.  we  esti¬ 
mated  the  unmodcled  water  vapor  error  as  approximately  6 
cm  rms.  This  is  a  global  estimate.  The  error  may  well  be 
somewhat  higher  in  the  Gulf  Stream  region.  We  ascribe  an 
estimate  of  6-10  cm  rms  for  the  total  of  both  higher-order 
water  vapor  and  ionospheric  range  error.  However,  the 
special  sensor  microwave  imager  (SSM/I)  based  studies  of 
Phoebus  and  Hawkins  [  1990).  while  carried  out  in  a  different 
region,  indicate  that  individual  instantaneous  passes  of  Geo¬ 
sat  ERM  altimetry  (such  as  the  simultaneous  pass)  can  be 
contaminated  with  even  larger  episodic  water  vapor  signa¬ 
tures. 

Nonlinear  errors  in  the  Schwiderski  tide  model  over  the 
relatively  short  track  segments  used  in  this  study  are  likely 
to  be  negligible  \Schwiderski.  1980).  while  again,  any  error 
possessing  a  linear  trend  alongtrack  will  be  implicitly  re¬ 
moved  during  the  TB  adjustment  (step  2). 

We  estimate  the  sum  of  the  error  due  to  surface  wave 
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TABLE  4.  Standard  Deviations  in  Alongtrack  Geoid  Profiles 


Standard  Deviation,  cm 


Ground  Track 

July  — 
April 

AXBT- 

Lamont 

AXBT- 

MM 

MM- 

Lamont 

A239 

27.76 

17.22 

83.74 

.4240 

15.69 

24.53 

36.72 

92.96 

31.00 

30.40 

A24I 

13.39 

37.27 

43.68 

122.1 

32.00 

47.19 

A242 

46.53 

34.46 

164.5 

Al 

81.39 

91.01 

258.7 

A2 

55.54 

85.93 

232.1 

A3 

16.24 

49.73 

97.74 

169.0 

53.83 

92.87 

A4 

69.95 

32.30 

123.8 

A5 

109.9 

62.32 

132.4 

Mean 

15.11 

51.62 

55.99 

153.26 

Differences 


MM.  Marsh/Mader  geoid  heights. 


skewness  and  residual  EM  bias  as  approximately  2  cm  [see 
Lybanon  and  Croat,  I987|. 

Finally,  as  will  be  seen  later  in  section  4.  residual  orbit 
determination  error  q  contaminates  the  estimated  geoid 
profiles  with  an  alongtrack  quadratic  with  an  amplitude  of 
approximately  8  cm. 

Regression  error  f(D,)  for  either  D15  or  DI7  is  estimated 
to  be  somewhat  larger  than  ±7  cm.  Error  in  representing  the 
total  topography  by  its  dynamic  (baroclinic)  component  is 
chiefly  due  to  the  presence  of  a  barotropic  component  in  the 
simultaneous  topography  ( 4 ,B I,, ) .  Using  IES/PG  data  from 
the  initial  (June  1985  to  July  1986)  REX  deployment.  Hallock 
t’t  ai.  [  1 989[  estimate  the  rms  barotropic  topographic  vari¬ 
ability  as  ±  10  cm.  We  accept  this  as  the  level  of  uncertainty 
due  to  the  simultaneous  barotropic  surface  topography 

Table  4  presents  standard  deviations  (i.e..  mean  differ¬ 
ences  removed)  between  various  geoid  height  estimates 
including  coincident  transects  through  the  Marsh/Mader 
\Mtir.sh  and  Chang.  1978]  and  Lamont  [  Wcsscl  and  Watts. 
1 988 1  gravimetric  geoids.  Refer  to  Figures  4  and  5  for  the 
location  of  transects. 

Independent  AXBT  sections  during  both  April  1987  and 
July  1987  were  collected  along  three  of  these  tracks  (A240. 
A24I.  and  A3;  see  Figure  4a).  Note  that  the  resulting  two 
estimated  geoid  profiles  are  independent  except  that  the 
same  ensemble  mean  sea  level  //  was  used  in  the  adjustment 
of  the  two  simultaneous  revs  for  orbit  error.  Ignoring,  for  the 
moment,  seasonal  steric  height  changes  from  April  to  July, 
differences  in  these  independently  estimated  geoid  profiles 
may  be  used  as  a  measure  of  the  total  error  of  our  procedure. 
The  average  standard  deviation  between  the  April  and  July 
estimates  for  these  three  ground  tracks  is  about  15  cm  (sec 
Table  4),  consistent  with  the  Hallock  cl  al.  1 1989)  assessment 
of  the  amplitude  of  barotropic  variability.  Biases  of  entire 
profiles  relative  to  each  other  appear  as  differences  between 
the  mean  sea  level  value  of  each  profile.  The  average 
profile-mean  difference  (July— April )  between  the  AXBT/ 
Geosat-estimated  geoid  profiles  (for  ground  tracks  A240. 
A24I,  anti  A3)  is  +  10.5  cm.  suggesting  that  seasonal  steric 
height  changes,  which  are  expected  to  be  of  approximately 


Fig.  4.  Ground  tracks  used  to  compute  geoid  profile  standard 
deviations  listed  in  Table  4  for  («)  repeats  of  A XBT  Geosal -derived 
geoid  profiles  and  (ft)  transects  through  the  Lamont  and  Marsh 
Mader  gravimetric  geoids  along  Geosat-F^RM  ground  tracks. 


this  magnitude  and  direction,  are  visible  in  our  estimated 
geoid  profiles. 

The  standard  deviation  (AXBT/Geosat-Lamont)  along 
tracks  shown  in  Figures  5«  and  5c  is  52  cm.  The  standard 
deviation  (AXBT/Geosat-Marsh/Mader)  along  tracks  shown 
in  Figures  5ft  and  5</  is  56  cm.  The  standard  deviation 
between  the  two  gravimetric  geoids  (Marsh/Mader-Lamont) 
along  tracks  shown  in  Figure  4ft  is  even  greater,  at  153  cm. 
Given  differences  in  reference  levels  between  the  two  gravi¬ 
metric  geoids  and  the  AXBT/Geosat-dcrived  geoid  profiles, 
as  well  as  the  arbitrary  adjustment  of  the  3000-dbar  level 
relative  to  the  altimetric  sea  level  in  our  estimation  of  the 
alongtrack  geoid  profiles,  it  is  meaningless  to  compare 
profile  mean  differences  between  the  various  geoid  profile 
estimates. 


4.  Gfosai  FRM  Orbit  Error  Anai  ysis 

To  compute  instantaneous  topography  C,  from  (4).  we 
would  like  to  justifiably  neglect  residual  short  wavelength 
orbit  error  in  the  reference  pass  (</,  ,„>sii<>ri  ■  The  residual 
orbit  error  in  the  NAG  ephenteris  is  very  long  wavelength 
(i.e..  wavelength  of  the  order  of  the  orbit's  circumference 
| see  Haines  ct  al..  I989|).  After  the  TB  adjustment,  the 
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Fig.  5.  Ground  (racks  used  (o  compute  geoid  profile  standard  deviations  listed  in  Table  4  for  (</)  differences  between 
April  19X7  AXBT'Geosat-derived  geoid  profiles  and  coincident  transects  through  the  l.amont  gravimetric  geoid.  I 
differences  between  April  19X7  AXBT'Geosat-derived  geoid  profiles  and  coincident  transects  through  the  Marsh 'Mader 
gravimetric  geoid.  (c  )  differences  between  July  19X7  AXBT'Geosat-derived  geoid  profiles  and  coincident  transects 
through  the  l.amont  gravimetric  geoid.  and  O/l  differences  between  July  19X7  AXBT  Geosat-derived  geoid  profiles  and 
coincident  transects  through  the  Marsh' Mader  gravimetric  geoid. 


residual  orbit  error  should  be  small  (less  than  10  cm:  see 
discussion  below).  The  extensive  number  of  alongtrack 
repeat  cycles  used  in  this  study  allow  s  for  an  assessment  of 
this  orbit  error.  The  adjustment  of  the  3000-dbar  level  to  an 
arbitrary  zero  relative  to  the  altimetric  reference  ellipsoid 
precludes  the  interesting  possibility  of  actually  calibrating 
the  orbit  error  by  using  the  AXBT  data. 

Since  the  long-wavelength  components  of  orbit  error 
dominate  fluctuations  in  the  cumulative  ensemble  mean  sett 
level,  we  may  assess  the  magnitude  of  this  error  by  exam¬ 
ining  the  stationarity  in  the  cumulative  value  of  //  from  ( 1 1  as 
N  increases.  As  an  increasing  number  of  repeated  revs  along 
a  given  ground  track  are  accumulated,  the  cumulative  rms  of 
the  biases  (f(/?)),  defined  here  as 

,v  ]  1 : 

X  m'-iN) 

•  i 


(w  here  B,  is  the  bias  of  rev  i.  i.e..  the  constant  term  from  the 
TB  adjustment,  and  V  is  the  total  number  of  revs  in  the 
ensemble)  approaches  its  stationary  value.  Again.  B,  repre¬ 
sents  the  mean  track  anomaly  from  tin  approximate  I -year 
average.  Table  5  lists  tlB)  for  each  of  the  13  ascending 
ground  tracks  in  the  REX  region.  The  cumulative  rms  of  the 
biases  ( f(/J))  may  be  used  as  an  estimate  of  the  amplitude  of 
the  long-wavelength  orbit  error  inherent  in  the  NAG  ephem- 
eris.  Our  a  priori  estimate  of  a  NAG  ephemeris  orbit  error  of 
approximately  3  m  rms  includes  both  geographically  corre¬ 
lated  and  uncorrelated  orbit  error.  However,  only  geograph¬ 
ically  correlated  orbit  error  is  visible  in  our  estimate  of  the 
alongtrack  topography.  Thus  as  can  be  seen  in  Table  3.  tin 
average  f(/J)  of  approximately  I  fib  cm  is  reasonable.  The 
long-wavelength  orbit  error  appears  at  shorter  wavelengths 
as  tin  artificial  alongtrack  gradient  in  sett  level.  This  gradient 
is  estimated  using  the  polynomial  expansion  of  a  weighted 
sum  of  a  sine  and  cosine  error  of  roughly  40.000  km 
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TABl.F.  5.  1  Ensemble  Orbil  luror  Statistics  for  N\V  Atlantic 

A  ml' ruling  Ground  Trucks 
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Track 

tiH  1. 

cm 
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cm 
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10 
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Ihh 

8.2 

10.3 

’Ground  track  for  which  there  is  no  calibration-quality  AXBT 
section. 
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Fig.  6.  Damping  as  a  function  of  number  of  repeal  revs  (sam¬ 
ples)  of  rms  mean  deviation  in  rev  biases  1  fi  in  section  4)  ( pluses  1. 
least  squares  lit  to  a  random  noise  curve  (squared,  and  estimated 
residual  short  wavelength  orbit  error  i/shori  (diamonds). 


wavelength  li  e.,  the  orbital  circumference): 

</  A ill  sin  l)  ■  C  cos  ID 

Air  •  Hit  (  ID  2!  Hti'  v.  •..  .1 

where  H  ZrrAt  I..  I.  is  orbital  circumference.  At  is  ground 
track  segment  length,  and  A  is  estimated  using  till). 

Table  5  lists  the  estimated  residual  orbit  error  after  re¬ 
moval  of  a  till  and  bias  li.e..  r/sttoKt  t,Z  -■*  over  ground 
track  segments  (At)  of  approximately  20(H)  km.  Note  that 
after  removal  of  a  tilt  and  bias,  sea  level  measurements  may 
still  he  contaminated  hv  approximately  8  cm  of  residual 
quadratic  orbit  error  over  these  segments.  Of  course,  the 
temptation  to  remove  increasingly  higher  order  polynomials 
from  sea  level  is  abated  by  the  knowledge  that  to  do  so  will 
increasingly  run  the  risk  of  inadvertent  removal  of  actual 
oceanographic  signal. 

As  the  geographically  uncorrelated  component  of  orbit 
error  is  a  once  per  rev  error,  it  should  be  decorrelated  over 
the  244-rev  interval  between  consecutive  repeats  of  a  given 
ground  track.  Thus  as  the  number  of  repeats  of  a  given 
ground  track  increases,  the  cumulative  rms  bias  or  mean 
deviation  ( fi).  given  by 

fi  1 1  r  l/f,l  rl II,  \  )  )  ""  I  /  -  I ) | 1  : 

should  monotonically  decrease  like  random  noise  as  I  (,\ 

I)  .  Note  that  as  defined  above,  the  rms  mean  deviation 
approaches  zero  at  /  -  ,Y.  For  example.  Figure  h  illustrates 
for  ground  track  A.3  the  time  history  of  (he  decrease  in  fi. 
slight  decrease  in  (directly  due  to  the  decrease  in  fi). 

and  a  least  squares  tit  of  a  FT,V  I)1-  or  random  noise 
curve  to  fi.  Similar  plots  for  each  of  the  1 .3  ascending  ground 
tracks  in  the  RF.X  region  are  characterized  by  a  fi  with  a 
sharper  rate  of  decrease  than  (he  least  squares  (it  of  a 
random  noise  curve. 

liquation  ( 1 1  is  used  to  establish  a  stationary  mean  sett 
level  value  within  each  bin  (//).  This  accomplishes  the 
damping  of  geographically  uncorrelated  orbit  error.  In  addi¬ 
tion.  other  time  variable  components  of  the  measured  sea 
level  (e.g..  variable  surface  topography  or  even  unmodeled 
variable  environmental  corrections)  tire  damped  in  the  cu¬ 


mulative  estimate  of  //.  Any  temporally  resolved  oceano¬ 
graphic  component  of  sea  level  (e.g..  surface  topography 
variability  of  frequencies  lower  than  the  sampling  frequency 
of  ae  collinear  altimetry  )  will  tend  to  damp  faster  than 
purely  random  noise  (see  appendix).  In  the  region  of  the  Gulf 
Stream,  surface  topographic  temporal  variance  can  be  as 
large  as  several  tens  of  centimeters  < i . e . .  only  somewhat 
smaller  than  the  magnitude  of  the  total  geographically  un¬ 
correlated  orbit  error).  As  is  shown  in  the  appendix,  collin¬ 
ear  repeats  which  sample  faster  than  the  significant  fre¬ 
quency  of  these  surface  topographic  fluctuations,  which 
chiefly  arise  from  the  propagation  of  Gulf  Stream  meanders, 
result  in  a  cumulative  ensemble  mean  sea  level  II  that  is 
uncontaminated  by  these  fluctuations.  Thus  the  rapid  de¬ 
crease  in  fi  suggests  temporal  resolution  of  a  significant 
portion  of  the  surface  topographic  variability  by  the  repeat 
altimetry.  Additionally,  the  time-incoherent  components  of 
surface  topographic  fluctuations  will  tend  to  result  in  the 
rapid  stationaritv  of  the  cumulative  ensemble  mean  II.  For 
some  ground  tracks  the  decrease  in  fi  is  not  purely  mono- 
tonic  and  is  characterized  by  low-amplitude  periodicity  after 
longer  (several  months)  averaging  periods,  suggesting  either 
aliased  or  low  frequency  (i.c..  time  scales  of  p  onths)  tem¬ 
porally  resolved  periodicities  in  sett  level.  See  Figure  7  for 
pertinent  examples  of  both  resolved  and  aliased  topographic 
fluctuations. 

An  important  parameter  associated  with  ensemble  aver¬ 
aging  of  collinear  sett  level  measurements  is  provided  by  the 
number  of  repeat  cy  cles  necessary  to  reduce  fl  to  tit  least  the 
level  iif  i/sit, >k i  ( i - e . .  the  number  of  repeat  cy  cles  necessary 
to  reach  stationaritv:  see  appendix).  Table  5  provides  this 
information.  Note  that  in  most  cases,  approximately  10 
repeal  cycles  (about  5-b  months)  must  be  accumulated 
before  the  ensemble  mean  orbit  error  is  reduced  to  the  noise 
floor  inherent  in  ti  I  B-only  adjustment  (i.c..  fi  is  reduced  to 
the  level  of  </s,lnK, ).  Beat  frequencies  between  the  collinear 
sampling  frequency  and  topographic  fluctuations  associated 
with  the  regional  oceanography  make  any  global  generaliza¬ 
tion  of  this  result  dangerous  and  lead  to  the  damping 
undulations  in  II  seen  m  Figure  7. 

The  ensemble  statistics  of  Table  5  and  the  associated 
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Kit;.  7.  Damping  within  an  averaging  hin  of  the  residual  rms  surface  topography  associated  with  the  propagation 
through  the  bin  of  a  simulated  100-em-amplitude.  330-km- wavelength.  8.4-cnvs  eastward  phase  speed  topographic 
fluctuation  (identified  by  Hiilliwcll  ami  .Wooers  119831  as  the  primary  mode  of  Gulf  Stream  meanders).  As  described  in 
the  appendix,  sampling  is  in  discrete  intervals  (At).  Damping  curves  for  repeat  periods  of  At  10.  70.  30.  40.  and  45 
days  tire  shown.  The  Geosat  KRM  has  a  repeat  period  of  17  days,  thus  resolving  the  major  mode  of  Gulf  Stream 
meander  fluctuations.  Severe  aliasing  of  this  45-day  meander  mode  occurs  for  repeat  periods  longer  than  77.68  days. 


curves  of  approach  to  stationary  (e.g..  Figure  6)  contain 
only  geographically  uncorrelated  orbit  error  (by  definition). 
However,  geographically  correlated  orbit  error,  which  is 
unaffected  by  the  averaging  process  of  (II.  is  a  potentially 
major  obstacle  to  the  realization  of  two-dimensional  maps  of 
surface  topography.  Such  maps  depend  upon  the  absence  of 
corrugations  or  biases  from  one  ground  track  to  the  next, 
regardless  of  w  hether  these  biases  arise  from  geographically 
correlated  or  uncorrelated  orbit  error. 

It  is  interesting  to  follow  the  propagation  of  both  geo¬ 
graphically  correlated  and  uncorrelated  orbit  error  through 
each  step  in  our  procedure  for  estimating  gcoid  profiles  from 
the  combination  of  AXBT  sections  and  simultaneous  Geosat 
ERM  altimetry.  A  summary  of  this  propagation  is  provided 
by  T  able  6.  where  the  algebraic  rules  of  error  propagation 


TABLE  6.  Slep-by-Step  Propagation  of  Geographically 
Correlated  and  Uncorrelated  Orbit  Error  in  the  AXBT  Geosat 
Geoid  Profile  Estimation  Procedure 


Step 

Description 

Correlated 

ok 

Contains 

Uncorrelated  OK 

1 

ensemble  mean  sea  level 

yes 

“no 

T 

TB  adjusted  instantaneous 
sea  level 

yes 

yes  (a  little)  no 

% 

AXBT  regression 
topography 

no 

ro 

4 

estimated  geoid  profile 

yes 

no 

surface  topography 
(computed  as  7— t) 

no 

no 

Rules  of  error  propagation  are  as  follows:  yes  yes  no:  yes  • 
no  no  •  yes  yes:  no  •  no  no;  ves  >  yes 


are  noted.  In  summary,  note  that  the  ensemble  mean  sea 
level  (step  1  in  section  2)  contains  only  geographically 
cot  related  orbit  error  (by  definition).  Adjustment  of  each 
instantaneous  sea  level  pass,  based  upon  removal  of  the 
linear  trend  (tilt  and  bias)  in  the  instantaneous  sea  level 
residual  with  respect  to  the  ensemble  mean,  propagates 
correlated  orbit  error  w  hile  introducing  a  very  small  amount 
of  uncorrelated.  higher-order  orbit  error  (due  to  quadratic 
and  higher-order  orbit  error  not  removed  from  the  instanta¬ 
neous  pass).  We  assume  that  the  dynamic  topography  based 
upon  the  AXBT  sections  contains  little  or  no  geographically 
correlated  error.  Hence  the  resulting  estimated  geoid  profile, 
w  hile  contaminated  by  geographically  correlated  orbit  error, 
is  not  subject  to  geographically  uncorrelated  orbit  error. 

Finally,  since  both  the  instantaneous  adjusted  pass  (pro¬ 
duced  in  step  2  of  section  2)  and  the  estimated  geoid  profile 
(produced  in  step  4  of  section  2)  contain  geographically 
correlated  orbit  error,  their  difference  (i.e. .  the  surface 
topography)  is  contaminated  by  neither  geographically  cor¬ 
related  nor  uncorrelated  orbit  error.  This  is  an  important 
characteristic  of  total  alongtrack  surface  topography  when 
computed  using  the  AXBT  Geosat-estimated  geoid  profiles. 
In  contrast,  absolute  surface  topography  computed  using  a 
gravimetrically  derived  geoid  profile  is  contaminated  by 
geographically  correlated  orbit  error.  While  this  contamina¬ 
tion  may  not  be  particularly  important  for  the  computation  of 
such  quantities  as  alongtrack  topographic  variability,  geo¬ 
graphically  correlated  orbit  error  can  introduce  large  track- 
to-track  corrugations  (in  the  cross-track  direction)  in  the 
inferred  surface  topography.  Unless  the  track-to-track  bi¬ 
ases  are  adequately  removed,  such  corrugations  would 
greatly  limit  (if  not  completely  preclude'  the  use  of  altimeter 
dala  for  objective  input  into  regional  numerical  models  of  the 
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ocean  circulation  (one  of  the  major  REX  objectives). 
Clearly,  care  must  be  taken  in  any  altimetric  attempt  to 
estimate  surface  topography  using  a  geoid  which  was  con¬ 
structed  independently  of  that  altimetry. 

5.  Conclusions 

Our  analysis  of  Geosat  ERM  collinear  altimetry  and 
simultaneous  AXBT  sections  has  resulted  in  the  following 
conclusions. 

1.  Estimated  alongtrack  geoid  profiles  across  the  Gulf 
Stream  appear  to  be  repeatable  to  within  10-20  cm  rms.  This 
repeatability  is  consistent  with  an  anticipated  minimum  error 
due  to  the  presence  of  any  instantaneous  barotropic  compo¬ 
nent  in  the  surface  topography  of  ±10  cm.  AXBT/Geosat- 
estimated  geoid  profiles  typically  agree  better  with  available 
gravimetric  geoids  than  the  several  available  gravimetric 
geoids  agree  with  each  other. 

2.  The  Geosat-ERM  NAG  ephemeris  radial  orbit  error 
over  the  NW  Atlantic  has  been  assessed  from  the  I -year  time 
series  of  repeated  measurements  of  sea  level.  The  rms  magni¬ 
tude  of  the  geographically  uncorrelated  orbit  error  of  approxi¬ 
mately  ±  1-2  m  is  in  agreement  with  earlier  simulation  studies 
designed  to  assess  errors  in  the  NAG  ephemeris.  Additionally, 
observed  orbit  error  scales  are  consistent  with  the  notion  that 
this  error  occurs  chiefly  at  a  frequency  of  once  per  revolution. 

A  more  desirable  variation  on  this  study  would  use  deep 
air  expendable  C'TDs  ( AXCTDs)  in  place  of  the  AXBTs.  The 
availability  of  reliable  deep  AXCTDs  would  allow  geoid 
profile  estimation  in  those  regions  where  the  dynamics  are 
not  clearly  dominated  by  temperature  (e.g..  NE  Pacific- 
subarctic  frontal  region  or  the  Greenland-lceland-Norwegian 
Sea  region).  The  marriage  of  satellite  altimetry  and  comple¬ 
mentary  air-dropped  expendable  instrumentation  anticipates 
a  long  and  fruitful  future. 

Appendix:  Defining  Mean  Topography 
With  Coi  i  inear  Ground  Trac  ks 

In  this  appendix  we  examine  the  rate  of  approach  to 
stationarity  of  the  cumulative  ensemble  mean  of  a  discretely 
sampled  fluctuation.  The  unavoidably  asynoptic  sampling  of 
surface  topographic  fluctuations  by  collinear  altimetry  will 
impact  the  nature  of  the  approach  of  the  ensemble  mean 
topography  to  stationarity  (e.g..  Figure  6).  In  general,  col¬ 
linear  repeat  periods  shorter  than  the  significant  time  scale  of 
the  topographic  fluctuations  will  resolve  these  fluctuations, 
and  the  resulting  approach  to  stationarity  by  the  ensemble 
mean  will  be  accelerated.  Alternately,  collinear  repeat  peri¬ 
ods  longer  than  the  significant  time  scale  of  the  topographic- 
fluctuations  will  alias  the  unresolved  fluctuations  and  result 
in  a  topographic  component  of  the  ensemble  mean  sea  level 
profile  which  does  not  monotonically  approach  stationarity. 
This  appendix  provides  detailed  discussion  of  this  issue 
using  as  an  example  the  cross-ground  track  propagation  of  a 
typical  Gulf  Stream  meander  across  a  fixed  collinear  Geosat 
ERM  groundtrack. 

Suppose  that  the  observed  sea  surface  topography  within 
an  averaging  bin  at  time  t  (give  by  £(/)>  is  equal  to  a 
fluctuating  component  plus  a  time  invaiiant  component: 

on  -  i"t.».  'F.  n  *  c„  (AD 

We  consider  the  time  variable  component  £'  at  position  v 
(the  cross-track  direction)  to  arise  due  to  a  topographic  wave 


(i.e.,  meander)  of  phase  'l/.  The  time  invariant  component  £0 
represents  the  true  temporal  mean  surface  topography.  We 
now  examine  the  rate  at  which  the  running  ensemble  average 
estimate  of  the  mean  topography  Z  approaches 
Without  any  loss  of  generality  we  may  consider  the  time 
variable  component  to  be  given  by  a  series  of  sine  waves,  so 
that 


M 

C’U,  'F.  /)  =  2  A„,  sin  [A,„(.v  -  ,v0 -<  „,/)  + 'F,„]  (A2) 

m  =  1 

where  we  consider  M  wave  components  each  of  wave 
number  km.  phase  speed  cm,  and  phase  Note  that  ,v 
represents  the  cross-track  direction,  where  ,v„  is  some  cross¬ 
track  reference  position  defining  the  location  of  the  ground 
track  bin.  The  “exact"  repeat  Geosat  ERM  ground  tracks 
are.  in  fact,  allowed  to  wander  over  a  cross-track  envelope 
of  up  to  ±1  km.  Thus  ,v  is  not  exactly  equal  to  .v„.  We 
consider  only  that  projection  of  wave  components  propagat¬ 
ing  transversely  across  track. 

Flying  in  an  exact  repeat  orbit,  the  altimeter  samples  the 
topography  at  point  x  in  some  regular  time  interval  given  by 
A t.  where  A t  =  17  “days"  for  a  17-day  repeat  orbit  (where 
a  “day”  equals  24  hours  +  the  daily  total  precession  rate  of 
the  satellite's  orbit).  Thus  in  the  case  of  the  Geosat  ERM 
17-day  orbit.  At  =  17  x  86651  s.  From  the  collinear  altimetry 
the  time  invariant  component  of  the  topography  £0  at  point  ,v 
on  a  collinear  ground  track  is  estimated  by  the  ensemble 
temporal  mean  Z  of  the  regularly  sampled  topography  at 
point  v.  given  by 

Z  =  f’(.v.  M'.  r)  +  fo  (A3) 

where  the  overbar  represents  a  running  ensemble  average. 
We  are  interested  in  the  time  rate  at  which  the  running 
ensemble  mean  of  the  fluctuating  component  £'  approaches 
zero  (i.e..  the  rate  at  which  Z  approaches  f„).  For  discrete 
sampling  at  intervals  of  A/,  this  is  given  by  the  series 

-  1  A 

t'l.v.  M'.  Nil)  =  -  2,  C(.v.  V.  nit) 

N  +  I  ‘ 


I 

N  +  I 


u 

zL  ‘3 m  sin  (Am(.\  —  .v0 

m  i 


-  c„,(. A.')  •  'F,„)|  (A4n) 

where  N  is  the  total  number  of  samples  in  the  ensemble  of 
tracks. 


We  desire  an  analytic  expression  for  £'  for  the  record 
length  Nit.  First  wx’  expand  the  expression  for  ('  as 


n.v.  'F.  Nit) 


I 


v 

V 


N 


V 


A„,  sin  (A,„(.v  -  .v (}) 


F  'I',,,)  cos  (A„,<  „,uAM| 


A  V 

—  ^  ^ .  A nl  cos  (A,rl(.\  -  Aid 

n  0  m  l 


'1^ fr,)  sin  (Aff,(  ninSt) 

( A5) 
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Recall  the  series  formulae 


^  sin  na  =  2  sir 


sin  [3  ( N  +  I )a]  sin  (a  S' a  ) 
sin  (a/2) 


v  sin  [( N  +  sla]  ^ 

X  cos  na  = - e - 

..  _ 2  sin  (a/2)  2 


Thus  (A5)  may  be  rewritten  as 


C'(x,  Silt) 


2  A„,  sin  [A ,„( ,v  -  ,v«)  +  'F„,]Fj(N.  k,„c„,St) 


X  A„,  cos  [A,,j.v  -  .v,,)  +  A„,<  „,Ar)  f  (A6) 


where  the  functions  F,  and  Fs  (for  M  orthogonal  wave 
components)  are  defined  as 

sin  [(A  +  i)A  A/]  , 

EjW,  kmcm\t )  =  — — — - - —  +  - 

2  sin  (kmc,„St/2)  2 


Fs(Af.  A ,„c„, A/)  — 


sin  [3  (A  +  l)A,„cmA/j  sin  (I  A/A„,<„,Ar) 
sin  (A,„r,„A//2) 


Wc  now  consider  only  one  of  the  M  components.  Addition¬ 
ally,  we  shall  at  this  point  neglect  the  effects  of  very  slight 
offsets  in  the  ground  track  location  and  simply  assume  that  v 
=  .v„  (i.e..  our  analysis  becomes  strictly  Eulerian).  We  are 
interested  in  the  rate  at  which  the  power  contained  in  the 
residual  £'  due  to  this  monochromatic  fluctuation  ap¬ 
proaches  zero  as  a  function  of  N.  In  general,  the  phase  V  is 
arbitrary,  since  we  will  lack  any  independent  realization  of 
the  topographic  wave.  Hence,  we  desire  a  phase-averaged 
rms  measure  of  this  tcsidual  power,  which  will  be  given  by 


<<T:>  =  — — ; -  {/V<VA'>  -  Fd/VA/)  +  F,(AA/)[  (A7) 

(/V  +  1)- 


— 1  *  f-’T  —  1 

<<r->  *  —  I  u"('i'.  aa/)]-</'P 

2w  Jo 

I  r:  rr  ,  I 

/»,( A 'Sn  -  —  I  F\A 2  sin  -  'I'  </'!'  -  - 

Jo 


I  Civ 

Pi(N±t)  s  —  I  h  xt\A  -  sin  1 

2jr  Jo 

FdA/An  -  —  I  /-;A‘  cos- 

2jr  Jo 


•■A  -  sin  'I'  cos  M'  (/'P  -  0 


Thus  (A7(  simply  becomes 


The  square  root  of  the  residual  power  (i.e..  residual  surface 
topographic  rms)  defined  by  (A8)  is  depicted  in  Figure  7  for 
various  altimeter  repeat  periods  when  sampling  the  most 
energetic  scale  of  Gulf  Stream  meanders  [see  Halliwell  and 
Mooers.  1983], 

Finally,  it  is  impoi  lain  to  note  that  this  fixed  point  (.v  -  .v„) 
analysis  made  no  attempt  to  take  advantage  of  cross-track 
information  on  the  phase  of  the  fluctuation.  In  actual  prac¬ 
tice.  if  the  total  mission  duration  is  longer  than  the  time 
required  for  the  topographic  fluctuation  to  propagate  from 
one  track  to  the  next  adjacent  track  (which  is  certainly  the 
case  for  Gulf  Stream  meanders  as  sampled  by  the  17-day 
repeat  Geosat  ERM  orbit),  then  cross-track  correlations  in 
the  topographic  time  series  might  be  used  to  obtain  a  fit  value 
for  A.  c.  and  VF.  Spatial  smoothing  of  the  data  will,  to  some 
extent,  accomplish  the  same  thing,  although  not  nearly  as 
precisely  or  eloquently. 
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